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(57) Abstract: A demodulator [500] for use in an ultrasound information processing system [200] is described, where in a preferred 
embodiment the demodulator [500] uses a half -band filter [508] to perform mirror-cancellation during the demodulation process, 
The demodulator comprises a mirror canceling quadrature mixer that mixes the input signals with half-band filters [508, 510] hav- 
ing a cutoff frequency of one-fourth the sampling fi^uency. In an ultrasound demodulator in wideband input signal is achieved. 
Additionally, the number and complexity of mirror canceling filter operations significantly reduced, allowing for easier hardware 
implementation using less expensive, off-the-shelf components. Furthermore, simplified hardware for digitally providing harmonic 
imaging applications is achieved. 
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DEMODULATING 
WTOE-BAND ULTRASOUND SIGNALS 

5 FIELD 

This patent specification relates to the field of ultrasound information processing 
systems. In particular, it relates to a method and system for the demodulation of 
ultrasound signals for use in ultrasound systems. 

10 BACKGROUND 

In recent decades ultrasonic imaging technology has played an increasing role in 
examining the internal structure of living organisms. The technology has applications in 
diagnosis of various medical ailments where it is useful to examine structural details in soft 
tissues within the body. Ultrasound imaging systems are advantageous for use in medical 

15 diagnosis as they are non-invasive, easy to use, and do not subject patients to the dangers 
of electromagnetic radiation. Instead of electromagnetic radiation, an ultrasound imaging 
system transmits sound waves of very high frequency (e.g., 2 MHz to 10 MHz) into the 
patient and processes echoes reflected from structures in the patient's body to form two 
dimensional or three dimensional images. 

20 More recently, ultrasonic imaging systems have been able to detect blood flow as 

well as tissue amplitude along the axis of the interrogating ultrasonic wave. These blood 
velocities are related by the Doppler effect movement within the body and blood flow is 
information of high diagnostic significance for certain diseases. This and other ultrasoxmd 
background information is disclosed in Zagzebski, Essentials of Ultrasound Physics 

25 (Mosby 1996), the contents of which are hereby incorporated by reference into the present 
disclosure. 

Signal demodulation represents a key preliminary step in converting reflected 
ultrasound signals into a usable representation on an output display. FIGS. 1 A and IB 
show exemplary plots of ultrasonic waves as transmitted into the body and received from 
30 the body, respectively, by an ultrasoimd transducer (not shown). In a typical ultrasound 
system operating in pulse-echo mode, the transducer transmits a signal 1 00 comprising a 
series of time-Iimited modulated bursts as shown in FIG. 1 A, and then receives a reflected 
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signal 102 as shown in FIG. IB. The plot 100 of FIG. lA is of a generic voltage signal 
which, as known in the art, is converted into sound waves by piezoelectric devices in the 
transducer for transmission into the body. The plot 102 of FIG. IB is likewise a generic 
voltage signal resulting from the conversion of reflected sound waves into voltages by the 
5 piezoelectric devices in the probe. 

It is to be appreciated for purposes of the present disclosure that system specifics 
such as frequency ranges, pulse durations, and the like are given by way of example only 
and are not intended to limit the scope of the preferred embodiments. By way of example 
and not by way of limitation, typical parameters for real-time ultrasound imaging 

10 applications would include: a carrier frequency Fc of 8 MHz (depending on the application, 
ultrasound implementations may use Fc selected from the range 2 MHz - 1 5 MHz); a burst 
repetition frequency BRF of 2000 Hz, Le„ 1 burst is sent every 500 |is (real-tune imaging 
systems may use a BRF selected from the range 2000 Hz - 4000 Hz); and a burst duration 
BD of 1 |is. As shown in FIG. IB, the reflected signal x(t) in plot 102 can be expressed as 

15 a sinusoid at the carrier frequency Fc modulated by an envelope signal A(t) and a time- 
varying phase term (|)(t), as can be expressed by Equations (l)-(3): 



x(t) = A(t)cos{a)ct + <|)(t) } { 1 } 

20 o)c = 27iFc {2} 

x(t) = Re{A(t)e'*^'>e^^"'} {3} 

Depending on the ultrasound implementation, it is the content of the envelope 

25 signal A(t) and/or the content of the phase signal (|>(t) that is of interest in generatmg the 

ultrasound system display output. For example, in B-mode imaging it is the envelope 

signal A(t) that is of interest because A(t) is proportional to the reflective index of the 

target tissue. As another example, in Doppler mode imaging it is the phase signal (|)(t) that 

is of interest because as the velocity of fluid flow in the target area is proportional to its 

30 first derivative d<|)/dt. Because die basic analytic signal A(t)e**^^^ modulated by the carrier 

frequency Fc in Equations (l)-(3) is usually found to contain frequency content in the 

-2- 



wo 01/54583 



PCT/DSOO/35723 



several-megahertz range, which is not significantly less than the carrier frequency Fc, the 
resultant signal x(t) is generally classified as a wideband signal. 

FIG. IC shows frequency spectra 104, 106, and 108 of the wideband signal x(t) 
corresponding to dififering exemplary target depths of 1 cm, 10 cm, and 20 cm, 
5 respectively. The spectra 104,106, and 108 correspond to the amplitude of the Fourier 
transform X(f) of the signal x(t) taken over small intervals of time around points 
corresponding to the expected arrival time of burst reflections from depths of 1 cm, 10 cm, 
and 20 cm, respectively. As illustrated by FIG. IC, the frequency spectrum X(f) is time 
varying in its characteristics. As known in the art, this is due to a frequency and depth 

10 attenuation factor of about 0.6 dB/MHz-cm for a typical ultrasound system on a typical 
human target, which causes the spectrum to be increasingly skewed toward lower 
frequencies in the far field (i.e., at greater depths) than in the near field {i.e., at lesser 
depths). Accordingly, there are different half-power bandwidths W of the analytic signal 
A(t)d*^^^^ as a function of depth, and in general Wdepth^di > Wdepth=d2 for d2>dl . For 

15 purposes of clarity of disclosure, and not for purposes of limiting the scope of the preferred 
embodiments, the analytic signal A(t)e^*^^^ is characterized herein as having a generic upper 
frequency limit W/2, where W is a half-power bandwidth for a depth in the near field. It is 
to be understood, however, that in practical systems a modified value for the upper 
frequency limit may be used without departing from the scope of the preferred 

20 embodiments. 

FIG. 2 shows a block diagram of a conventional ultrasoimd system 200 in 
accordance with prior art. Ultrasound system 200 comprises a transducer 202, a front end 
processor 204, a demodulator 206, and an amplitude detector and display apparatus 208. 
As described supra, transducer 202 transmits focused acoustic signals into the body and 

25 sends an analog response signal to front end processor 204. Front end processor 204 
performs preliminary operations such as depth gain compensation and then transmits the 
result, commonly called the RF signal, to demodulator 206. For purposes of clarity and 
simplicity, the analog signal x(t), shown as plot 102 of FIG. IB, is referred to herein as the 
RF signal, with depth gain compensation and/or other preliminary processing having 

30 already been performed by front end processor 204. 
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In the system of FIG. 2, the demodulator 206 perfonns processing steps directed to 
extracting the envelope signal A(t) and the phase signal ^(t) from the RF signal x(t), and 
then transmits the result to amplitude detector and display apparatus 208 for downstream 
processing and eventual image display. In a simplest form, for B-mode processing such a 
5 demodulator could simply comprise an analog envelope detector having a full-wave 
rectifier and a low-pass filter, as described in Carlson, Communication Systems: An 
Introduction to Signals and Noise in Electrical Communication (McGraw-Hill, 3^** ed. 
1986), the contents of which are hereby incorporated by reference into the present 
disclosure. However, in most practical ultrasound systems in which both the amplitude 

10 signal A(t) will be desired for a B-mode processing mode, and where phase signal <t)(t) will 
be desired for Doppler processing mode, a coherent detection scheme is used wherein 
demodulator 206 comprises a quadrature mixer. While older prior art implementations 
were analog in nature, most newer implementations are digital. For purposes of the 
present disclosure, and without loss of generality, it may be presumed that front end 

15 processor 204 incorporates an anti-aliasing filter followed by an analog to digital converter 
operating at a sampling frequency of Fj. 

FIG, 3 shows a block diagram of demodulator 206 in accordance with prior art. 
Demodulator 206 comprises a local oscillator 302, mixers 304 and 306, and low pass 
filters 308 and 310. In accordance with conventional digital signal processing principles, 

20 the sampling frequency Fs is at least satisfying the Nyquist rate, which for the system of 
FIG. 3, 2(Fc + W/2) is selected where, as described supra^ Fc is the ultrasound carrier 
frequency and W/2 is the upper frequency limit W/2 of the analytic signal A(t)e**^*^ Upon 
sampling, the digitized RF signal x(kT) results, with the sampling period being T = 1/ Fs, 
the signal x(kT) simply being cast herein as x(k). Local oscillator 302 is designed to 

25 generate sinusoids at a mixing frequency Fx that are in quadrature phase with each other, 
e.g., cos(27tFxk) and sin(27tFxk). Mixers 304 and 306 multiply the digitized RF signal x(k) 
with the respective quadrature-phase signals from the local oscillator 302, and the products 
are sent to low-pass filters 308 and 3 1 0, respectively. Low-pass filters 308 and 3 10 are 
designed to have sharp rolloffs at the upper frequency limit W/2 of the analytic signal 

30 ACt)^*^'^ 
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As known in the art, for Doppler processing the mixing frequency Fx is fixed at a 
frequency identical to the carrier frequency Fc. Using relationships known in the art and as 
described generally in Carlson, supra, in such case the demodulator 206 generates the 
signals I(k) and Q(k), where I(k) and Q(k) are related to the amplitude signal A(t) 
5 (expressed herein as A(k) in the digital domain) and phase signal <|)(t) (expressed herein as 
<t)(k) in the digital domain) as shown in Equations (4)-(8) below: 

A(k)e'*^^ = I(k)+jQ(k) {4} 

10 I(k) = A(k)cos<|»(k) {5} 

Q(k) = A(k)sin(Kk) {6} 

A(k) = V{l2(k) + Q2(k)} {7} 

(|)(k)=tan-'{Q(k)/I(k)} {8} 



15 



For B-mode processing, where it is desirable to obtain amplitude A(k) and where 
detection of phase <|>(k) is not important, the mixing frequency Fx is selected to be time- 

20 varying for increasing the signal- to- noise ratio. In particular, the mixing frequency Fx is 
selected to correspond to an instantaneous center frequency of the signal x(t) which, as 
shown in FIG. IC, becomes increasingly lower than the carrier frequency Fc as the field 
depth increases. In such case, the oscillator 302 is a swept frequency oscillator, with the 
mixer firequency Fx varying from a maximum value to a minimum value over each burst 

25 period BP. As known in the art, using a swept oscillator for the mixing frequency Fx 
provides better signal-to-noise performance and therefore better detection of A(k) as 
compared to using a mixing frequency Fx fixed at the carrier frequency Fc. However, it 
can also become impractical to recover the phase signal (|)(k) using this method. 

By way of nonlimiting example, and fiirther to the example of FIGS. IC, for B- 

30 mode processing the mixing frequency Fx may sweep from a maximum of 7.5 MHz in the 
near field to 5.0 MHz in the far field. It can be shown for such B-mode processing that the 
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amplitude signal A(k) can still bexomputed from the demodulator 206 outputs I(k) and 
Q(k) from Equation (7) above while, as discussed supra, the phase signal <|>(k) is generally 
not computed. 

FIG. 4, comprising parts 4-1 and 4-2. shows Fourier spectrum plots corresponding 
S to the operation of prior art demodulator 206 of FIG. 3. FIG. 4 shows a first plot 402 of 
the Fourier transform X(f) of the RF signal x(k), a second plot 404 of the Fourier transform 
of the local oscillator output cos(ci)xkT), a third plot 406 of the Fourier transform of the 
mixer 304 output, a fourth plot 410 of the Fourier transform of the local oscillator output 
sin(a)xkT) as multiplied by j, and a fifth plot 412 of the Fourier transform of the mixer 306 

10 output as multiplied by j. Note that Fourier transform plots 402, 406, and 412 repeat at 
intervals of Fs as they represent digital signals at a sampling frequency of Fs, and only the 
frequency range of -Fj to Fj is shown. 

Several problems, however, arise in ultrasound systems using the prior art 
wideband demodulator scheme of FIGS. 2 and 3. As shown on spectrum plots 406 and 

15 4 1 2, the low-pass filters 308 and 3 1 0 must each have a frequency characteristic similar to 
the filter characteristic 408 for appropriate mirror rejection, i.e. cancellation of the 
undesired mirror portion of the mixed signals that occur at the output of mixers 304 and 
306, respectively. As shown in FIG. 4, filter characteristic 408 of low-pass filters 308 and 
3 10 is required to have a very sharp roll-off for proper construction of the signals I(k) and 

20 Q(k), respectively. In FIG. 4, for example, it is shown that the distance between the 

desired signal and its undesired mirror signal is approximately (Fs-2Fc-W) if Fs<2(Fc-i-Fx), 
and approximately 2(Fx-W) if Fs>=2(Fc+Fx). Using the exemplary parameters of FIG. 4, 
this distance is only 32 - 2(8) - 6.4 = 9.6 Mhz. Implementing this sharp roll-off, however, 
can create large group delay distortion at the filter band edge that degrade images produced 

25 from these signals and cause ringing. Additionally, using prior art ultrasound sampling 
and filter frequencies, the computations for implementing the low pass filters 308 and 3 1 0 
are complex and time consuming, requiring general purpose digital signal processors of a 
very high order and very high power, Disadvantageously, the hardware needed to 
implement the prior art demodulator suffers from increased size, heat dissipation 

30 requirements, cost, and complexity due to the need for the sharp roll-off characteristic of 
low pass filters 308 and 310. 
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An additional problem with prior art demodulators as in FIG. 3 arises in the context 
of B-mode processing, which usually incorporates a swept mixing frequency Fx at local 
oscillator 304 as described supra. Because the cutoff frequency of each of the low pass 
filters 308 and 3 1 0 is fixed at a single value, the mirror cancellation performed by low pass 
5 fihers 308 and 3 10 is inherently better for certain values of mixing frequency Fx than 
others. A undesirable design choice is forced upon the ultrasound hardware designer. 
Given this dilemma, the design choice is usually made to select the cutoff frequency of the 
filters 308 and 310 such that near field reflections receive better mirror-canceling 
performance than far-field reflections. As a result of this design choice, there is inferior 

10 demodulation of far-field reflections as compared to near-field reflections. 

One attempted prior solution to the above problems is described in U.S. Pat. No. 
5,482,044, the contents of which are hereby incorporated by reference into the present 
application. The method described therein, however, which includes the use of a Hilbert 
transform, suffers from practical implementation difficulties and complexities, including 

15 problems related to the fact that the Hilbert transform is a noncausal signal which is not 
realizable in a physical system and can only be approximated. 

A further problem with the prior art ultrasound information processing system of 
FIGS. 2 and 3, as well as U.S. Pat. No. 5,482,044, supra, is found in the context of 
harmonic imaging. As known in the art, in harmonic imaging mode an ultrasound 

20 information processing system attempts to recover second harmonic information near 2Fc. 
The above cited prior art systems, however, provide no easy digital solution to the 
recovery of the second harmonic information, and in field implementations there is usually 
an analog high pass filter incorporated into the front end processor 204 that removes, in the 
analog domain, the information near Fc. The requirement to include an adjustable analog 

25 filter in the front end processor 204 with sufficient high pass characteristics presents 
problems of cost and complexity that it would be to desirable to avoid. 

Accordingly, it would be desirable to provide an ultrasound demodulator having an 
improved output with better mirror canceling, less group delay distortion, and less ringing. 
It would be fiirther desirable to provide an ultrasoimd demodulator having an 

30 improved output while also being easier to implement in hardware, using lower cost, oflF- 
the-shelf components with lesser footprint and heat dissipation requirements. 

-7- 
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It would be still further desirable to provide an ultrasound demodulator that, in 
amplitude-only implementations having swept-frequency mixing at a local oscillator, 
produces improved output in the far field as well as the near field. 

It would be still further desirable to provide an ultrasound information processing 
S system having the capability of performing harmonic imaging using low-cost, off-the-shelf 
components that implement a digital technique incorporated into the ultrasound 
demodulator hardware, and not requiring an analog high pass filter at a fi-ont end processor 
of the ultrasound information processing system. 

10 SUMMARY 

A demodulator for use in an ultrasound information processing system is provided, 
where in a preferred embodiment the demodulator uses a half-band filter to perform 
mirror-cancellation during the demodulation process. The demodulator can comprise a 
mirror canceling quadrature mixer that mixes the input signal with quadrature phase 

15 sinusoids oscillating at one-fourth of the sampling fi-equency, and then low-pass-filters the 
results with half-band filters having a cutoff frequency equal to one-fourth the sampling 
frequency. It has been found that in an ultrasound demodulator in accordance with a 
preferred embodiment, a significantly sharper mirror cancellation on the wideband input 
signal is achieved. Additionally, the number and complexity of mirror canceling filter 

20 operations is significantly reduced, allowing for easier hardware implementation using less 
expensive, off-the-shelf components. Furthermore, simplified hardware for digitally 
providing harmonic imaging applications is achieved. 

An ultrasound information processing system for which the demodulator is suitable 
can comprise an ultrasound transducer generating acoustic bursts at a carrier fi-equency, the 

25 ultrasound transducer receiving acoustic echo signals and generating an input signal 
therefi-om. The ultrasound information processing system filters the input signal with an 
anti-aliasing filter and digitizes the input signal at a sampling frequency to generate digital 
samples. A demodulator in accordance with the preferred embodiments comprises a 
quadrature mixer adapted to multiply the digital samples with quadrature-phase sinusoids 

30 at one quarter of said sampling frequency and to filter the results with half-band filters, 
each having a cutoff fi-equency of one-quarter of the sampling frequency, to produce 
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FIG. 6 comprises parts 6-1 through 6-4 and frequency spectra corresponding to the 
operation of the demodulator of FIG. 5; 

FIG. 7 is a block diagram of a demodulator in accordance with another preferred 
embodiment; and 

5 FIG. 8 illustrates frequency spectra corresponding to a harmonic imaging mode of 

the demodulator of FIG. 7 in accordance with a preferred embodiment. 

DETAILED DESCRIPTION 

FIG. 5 shows a block diagram of a demodulator 500 for use in an ultrasound 

10 information processing system in accordance with a preferred embodiment It is to be 
appreciated that in the following description, specific parameters that are presented are 
done so by way of example so as to help illustrate and describe preferred embodiments, 
and are not presented by way of limitation. Demodulator 500 optionally comprises in an 
anti-aliasing filter (not shown) and an analog-to-digital converter (not shown) operating at 

15 a sampling frequency Fj. The anti-aliasing filter and the analog-to-digital converter may 
altematively be considered as part of a front end processor located in front of the 
demodulator 500. Demodulator 500 comprises a first stage local oscillator 502, mixers 
504 and 506, half-band filters 508 and 5 10, a complex mixer 512, and a secondary local 
oscillator 514. Inputs to demodulator 500 include the wideband signal x(k) that has been 

20 sampled at a sampling frequency of Fs, and outputs include an in-phase component signal 
I(k) and a quadrature-phase component signal Q(k). The desired envelope signal A(k) is 
related to I(k) and Q(k) by equation (7) supra and, for Doppler mode processing, the phase 
signal (|)(k) is related to I(k) and Q(k) by equation (8) supra. The sampling frequency Fs is 
selected such that it is at least twice as great as the highest frequency in the wideband 

25 signal x(t). In a preferred embodiment, the sampling frequency Fj is 48 MHz. 

An anti-aliasing filter (not shown) is used to filter out unwanted noise or 
unnecessary signals from the input signal x(t) above the Nyquist frequency of ahalog-to- 
digital converter (not shown) in order to prevent aliasing in the analog-to-digital converter. 
The Nyquist frequency is one half the sampling frequency, i.e., The cutoff frequency 

30 of the anti-aliasing filter is preferably below 15 MHz but must be large enough to 

accommodate the information in the signal x(t). Using the exemplary parameters of Fs = 

-10- 
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48 MHz, Fc = 8 MHz, and W = 0.8Fc as described in this patent specification, a cutoff 
frequency of 1 5 MHz is allowable as the upper frequency limit of the signal x(t) is not 
above 1 1 .2 MHz. The output of the anti-aliasing filter is supplied to the analog to digital 
converter where it is sampled at the sampling frequency Fs, and the sampled signals are 
5 digitized. The digitized samples are denoted herein as x(k). 

As shown in FIG, 5, demodulator 500 comprises a first stage local oscillator 502 
for producing quadrature sinusoids at one-fourth the sampling frequency, ie. at Fs/4. 
According to a preferred embodiment, the frequency Fs/4 of the first stage local oscillator 
502 is advantageously selected such that mirror canceling operations may be subsequently 
10 performed by half band filters. The quadrature sinusoids at the output of the first stage 
local oscillator 502 are then multiplied by the signal x(k) at mixers 304 and 306, 
respectively. The formulae for the sinusoids generated by first stage local oscillator 502 
are given below in Equations. (9)-(I0): 

15 cos{27c(Fs/4)(kT)) = cos(k7i/2) = {1,0, -1, 0, 1, 0, -1, 0, . . .} {9} 

. sin{2K(Fs/4)(kT)) = sin(k7i/2) = {0, 1, 0, -1, 0, 1, 0, -1, . . .} {10} 

As reflected in equations (9)-(10), the use of Fs/4 as the oscillation frequency for 
20 the first stage local oscillator 502 results in an especially simple operation for the digital 
mixers 504 and 506 that avoids the need for any actual multiplications. Since the oscillator 
signal is zero for every other term of the cosine and sine components as shown in 
Equations (9)-(10), the output of each mixer 504 and 506 can automatically be decimated 
to one-half the sampling frequency, Fs/2, with a simple alignment in time delay. Further, 
25 since the non-zero values of the oscillator signal alternate between 1 and -1, the mixing of 
these values with the sampled signals can be accomplished simply by inverting every other 
sample. 

The signals from mixers 504 and 506 are fed to first stage half-band filters (HBF) 
508 and 510. Advantageously, in accordance with a preferred embodiment, the use of a 
30 half-band filter having a passband between -Fs/4 and Fs/4 allows for further simplified 
hardware. In particular, for digital finite-impulse-response (FIR) filter implementations. 
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the impulse response of a half-band filter is zero for every other impulse response sample. 
The resulting hardware efficiencies allow for extremely high performance FIR filters 
having very linear, sharp-rollofif characteristics to be implemented for half-band filters 508 
and 510 that are, at the same time, of a modest cost and complexity. In a preferred 
5 embodiment, a 67-tap FIR implementation is used having the remarkable rolloff 
characteristic that at 0.29Fs (only 0.04Fs above the cutoff frequency of 0.25 Fs), the 
attenuation is 90dB. 

It has been found that the use of a particular off-the shelf chip normally used in the 
telecommunications industry, and not formerly associated with usage in ultrasound 

10 systems, provides excellent results. The HSP 43216 chip is an off-the-shelf chip sold by 
Intersil Corp. of Melbourne, Florida, and has been found particularly adaptable for use an a 
demodulator according to the preferred embodiments. A single HSP 43216 chip can be 
used to implement the combination of the first stage local oscillator 502, the mixers 504 
and 506, and the half-band filters 508 and 510. The datasheet for the HSP 43216 chip, 

15 Intersil File No. 3365.7 (January 1 999), can be found on the Internet at 

www.intersil.com/data/fii/fii3/fii3365/fii3365.pdf, and application notes relevant to the 
HSP 43216 chip, Intersil File No. AN9063.2 (January 1999), can be found at 
www.intersil.com/data/an/an9/an9063/an9603.pdf The contents of these documents are 
hereby incorporated by reference into the present disclosure. 

20 As shown in FIG. 5, the outputs of half band filters 508 and 5 1 0 are denoted as l\(k) 

and Qi(k), respectively, as they represent intermediate signals in the generation of the 
demodulator 500 output. Demodulator 500 further comprises a complex mixer 512 and a 
second stage local oscillator 514. Upon the creation of the signals l\(k) and Qi(k) at the 
outputs of half-band filters 508 and 510, these signals are then fed into the complex mixer 

25 512 for rotation to the desired baseband signals I(k) and Q(k). The rotation is performed 
by mixing the signals Ii(k) and Qi(k) with quadrature phase sinusoids, shown as cosine and 
sine in FIG. 5, that are generated by secondary local oscillator 514 at an advantageously 
selected mixing frequency of Fs/4 - Fx- Here, Fx may be a constant or varying frequency 
fimction, and is generally the same mixing frequency that would be used in the single-stage 

30 demodulator of FIG. 2. In other words, for Doppler processing Fx is fixed at Fx = Fc, and 
for B-mode processing Fx is set to a swept frequency function, in particular an 

-12- 
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instantaneous center frequency of the wideband input signal x(t) for better signal-to-noise 
performance. 

Complex mixer 512 rotates the intermediate signals Ii(k) and Qi(k) to the desired 
baseband signals I(k) and Q(k) by performing operations that accomplish the complex 
5 multiplication shown in Equation (11): 

I(k) + jQ(k) = PiCk) + jQi(k)][cos{(G)s/4.G)x)kT} - jsin{(tDs/4-fi>x)kT}] {11} 

In particular, in performing operations that accomplish the complex multiplication 
10 above, complex mixer 512 generates the in-phase component I(k) by performing the 
multiplies and adds shown in Equation (12), and generates the quadrature component by 
performing the multiplies and adds shown in Equation (13): 



15 



I(k) = Ii(k)cos{(a)s/4-G)x)kT}+ Qi(k)sin{(a)s/4-a)x)kT} {12} 
Q(k) = Qi(k)cos{(c)s/4-a)x)kT} - Ii(k)sin{(Q)s/4-a)x)kT} {13} 



Digital hardware for implementing complex mixer 512 generally comprises 
multiplier-accumulator circuitry for performing the above multiplies and adds which, in 

20 accordance with a preferred embodiment, are performed on 16-bit values and generate 16- 
bit results. One suitable hardware chip is the LMAlOlO multiplier/accumulator chip 
available from Logic Devices. It is found to be advantageous for secondary local oscillator 
514 to be implemented as a lookup table comprising a static random access memory 
(SRAM) to generate the sine and cosine signals at a frequency of (Fs/4-Fx). In a preferred 

25 embodiment, the SRAM is 16K in size and stores 16 bit output values. 

In operation, the SRAM address counter within secondary local oscillator 5 14 is 
reset for each transmit synchronization burst and then incremented at a frequency of Fs 
during that burst interval. For Doppler mode operation, the SRAM simply provides the 
required sinusoidal values corresponding to a frequency of Fs/4 - Fx, where Fx is fixed at 

30 Fc- For B-mode imaging operations, the SRAM provides swept sinusoidal data for the 
depth frequency compensation ftinction, z.e., for creating the appropriate mixing frequency 
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Fs/4 - Fx where Fx sweeps between maximum and minimum values during each burst 
interval. 

Advantageously, in addition to providing Doppler and B-mode functionalities, the 
same SRAM and complex mixer hardware may be used (in conjunction with downstream 
5 low pass filter hardware, to be described injra) to achieve harmonic imaging mode 

functionality. This is accomplished by changing an offset frequency of the secondary local 
oscillator such that the frequency Fx is swept at values near 2Fc which, as will be described 
further infra^ allows information near the second harmonic of the carrier frequency to be 
isolated and observed. 

10 FIG. 6, comprising parts 6-1 through 6-4, shows spectrum plots corresponding to 

the operation of demodulator 500. FIG. 6 shows a first plot 602 of the Fourier transform of 
the RF signal x(k), a second plot 604 of the Fourier transform of the first stage local 
oscillator cosine term, and a third plot 606 of the Fourier transform of the mixer 504 
output. FIG. 6 further shows a fourth plot 610 of the Fourier transform of the first stage 

15 local oscillator sine term, and a fifth plot 612 of the Fourier transform of the mixer 506 
output. Upon being shifted in both directions by the intermediate frequency Fi, which is set 
to Fs/4 according to a preferred embodiment, the mixed signals 606 and 612 each contain a 
gap between the baseband portion and the undesired mirror portion approximately equal to 
(2Fc - W) if Fc<=Fs/4, and approximately equal to (Fj - 2Fc - W) if Fc>Fs/4. In an example 

20 in which Fc = 8 and W = 0.8Fc, this gap is about 9.6 MHz. Importantly, however, despite 
the fact that this gap is similar to the gap at the mixer outputs of FIG. 3 in a similar 
example described supra, there is superior mirror cancellation according to the preferred 
embodiments. This is at least in part because each gap is symmetric around Fs/4, and 
therefore a half-band filter can be used having an extremely sharp roUoff at Fs/4, The 

25 passbands of the half-band filters 508 and 510, showing an exemplary steep rolloff, are 
shown as elements 608 superimposed on the plots 606 and 612 of FIG. 6. 

Also shown in FIG. 6 are frequency spectrum plots corresponding to the operation 
of complex mixer 512 of demodulator 500. FIG. 6 shows a frequency plot 614 of the 
output of half band filter 508, which corresponds to the intermediate in-phase signal Ii(f), 

30 and a frequency plot 620 of the output of half band filter 5 1 0 multiplied by j, which 
corresponds to the intermediate quadrature phase signal Q\(f). Frequency plots 616 and 
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622 are also shown corresponding to the Fourier transforms of the outputs of secondary 
local oscillator 514. Additionally, frequency plots 618, 624, 628, and 630 are shown 
corresponding to the Fourier transforms of the intermediate multiplications performed in 
the complex mixer 512 indicated in Equations (12) and (13), supra. Finally, FIG. 6 shows 
S a plot 626 of the Fourier transform of the in phase component output I(k) of complex 
mixer 512, and a plot 632 of the Fourier transform of the quadrature phase component 
output Q(k) of complex mixer 512. Advantageously, as can be seen by comparing plots 
626 and 632 of FIG. 6 with plots 406 and 412 of FIG. 4, respectively, the frequency 
content of the signals at the baseband becomes identical to the desired baseband 

10 components for I(k) and Q(k) because of the selection of (Fs/4-Fx) as the frequency for 
secondary local oscillator 514. Additionally, because the mirror canceling was performed 
by very sharp half band filters, the signals 626 and 632 of FIG. 6 represent "cleaner" 
outputs than the signals 406 and 412 of FIG. 4, and are also achieved using less expensive 
and less complex hardvdrare in accordance with the preferred embodiments. 

15 FIG. 7 shows a demodulator 700 for use in an ultrasound information processing 

system in accordance with another preferred embodiment. Demodulator 700 comprises 
substantially the demodulator 500 of FIG. 5, to which is coupled second stage low pass 
filters 702 and 704, along with integrate-and-dump matched filters 706 and 708, 
respectively. The additional components of demodulator 700 provide for fiirther optimized 

20 output quality during Doppler and B-mode processing as described herein, and additionally 
provide for harmonic imaging mode processing when used in conjunction v^th properly 
offset frequencies at the secondary local oscillator 5 14 as described infra. As shown in 
FIG. 7, the outputs of demodulator 500 are relabeled as Ii(k) and Qi(k) to indicate that 
further processing is being performed to generate final demodulator outputs Io(k) and Qo(t) 

25 at the outputs of integrate-and-dump matched filters 706 and 708, 

In a preferred embodiment, dual programmable low pass filters are used to 
implement the low pass filters 702 and 704. In hardware implementations, it is found 
usefiil to use the HSP 43168, a programmable 8-tap FIR filter available from Intersil Corp., 
supra. Such programmable filters allow for the shifting of the cutoff frequency as needs 

30 arise. During Doppler mode and B-mode imaging, the low pass filters 702 and 704 are set 
to appropriate bandwidths to provide further noise reduction and attenuation of spurious 

-15- 



wo 01/54583 



PCT/USOO/35723 



high frequencies. The low pass filters 702 and 704 will have cutoff frequencies that 
depend on the frequency of the signal received. For example, when an increase in depth is 
required, the receiving mean frequency is shifted down. Correspondingly, the bandwidth of 
the low pass filters 702 and 704 should be reduced. 
5 FIG. 8 shows power spectra of the in phase signal component when propagating 

through the demodulator 700 in harmonic imaging mode in accordance with a preferred 
embodiment. Advantageously, the low pass filters 702 and 704, even though implemented 
using relatively inexpensive 8-tap FIR filters, can be readily adapted for use in conjunction 
with properly offset frequencies at the secondary local oscillator 514 to achieve harmonic 

10 imaging mode. Plot 802 of FIG. 8 shows the Fourier transform of a typical RF signal x(k) 
received during harmonic imaging mode, the RF signal x(k) having its first harmonic at a 
carrier frequency of 3.5 MHz (shown as plot element 804 in FIG. 8) and a second 
harmonic at 7 MHz (shown as plot element 806 and FIG. 8). In accordance with a 
preferred embodiment, the secondary local oscillator 514 generates frequencies of (Fs/4- 

15 Fx), either swept or fixed in nature, such that Fx is near the second harmonic 2Fc. Using 
relationships established with respect to the description of demodulator 500 supra, it can 
be shown that the plot 808 results for in phase component Ii(k) at the output of complex 
mixer 512. For simplicity and clarity of disclosure, only the in phase components are 
described with respect to harmonic imaging mode, it being understood that similar 

20 corresponding frequency contents will be experienced with die quadrature phase 
components. As shovm in plot 808, it is the second harmonic components that are 
congregated near the origin of the frequency plot, while the first harmonic components are 
piished out to +Fc and -Fc. 

In accordance wdth a preferred embodiment, for harmonic hnaging mode the low 

25 pass filters 702 and 704 are programmed to have the spectrum 810 identified as LPFh(f) in 
FIG. 8 having a passband at the origin and first notches at -Fc and +Fc. Advantageously, 
even though the low pass filters 702 and 704 are simple 8-tap FIR filters, the spectrum 810 
is readily achieved having distinct zero pass notches at — Fc and Fc. As shown in plot 812, 
which is a plot of the Fourier spectrum of signal hCk) at the output of low pass filter 702, 

30 the first harmonic components of the RF signal x(k) are effectively removed and the 
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second harmonic components are passed, thus achieving harmonic imaging mode 
functionality. 

As shown in FIG. 7, the signals tOc) and Qi(k) at the outputs of low pass filters 
702 and 704, respectively, are fed into integrate-and-dump matched filters MF 706 and 
S 708, respectively. Integrate-and-dump matched filters are known in the art, and a 

description can be found in Couch, Digital And Analog Communications Systems, 3rd ed. 
(MacMillan 1990) at pp. 497-503, the contents of which are hereby incorporated by 
reference into die present application. Generally speaking, the integrate-and-dump 
matched filters 706 and 708 provide for increased signal-to-noise ratio through matched 

10 filtering using an assumption that the transmitted pulse is a square wave for the burst 

duration BD (see FIG. 1 A). Additionally, the integrate-and-dump matched filters 706 and 
708 provide for an integration or averaging function of the data for subsequent image 
display, as demodulation operations at the sampling frequency Fs in accordance with the 
preferred embodiments results in more data than needed for visual display on most current 

15 ultrasound display equipment. Integrate-and-dump matched filters 706 and 708 are 
adjusted dependent on the display pixel density desired. For example, the filtered output 
may be further filtered at a lower cut-off frequency to prevent spatial aliasing on the 
display screen. With such an implementation, contrast resolution on the display screen is 
also improved. 

20 In accordance with another preferred embodiment, the demodulator signals I(k) and 

Q(k) at the output of demodulator 500 or 700 can be used by downstream circuitry in B- 
mode imaging to produce useful flow information, as follows. In particular, a phase 
detection circuit performs a tan'^ [Q(k)/I(k)] operation to determine phase information of 
the demodulated signal. Due to the large number of bits each signal has (16 bits each in a 

25 preferred embodiment), each signal is fu-st compressed Jfrom a 16 bit number to an 8 bit 
number by going through a RAM look up table to implement the tan"^ [Q(k)/I(k)] function. 
An amplitude detection circuit likewise computes the amplitude using Equation (7). The 
instantaneous phase and amplitude for each sample for an entire ultrasound frame are 
stored in memory as the next frame is being processed. The phase information at each 

30 sample position is then correlated (or compared) with the same sample position in the last 
frame. If there is blood flow at that sample position, the phase will be different from frame 
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to frame. The stationary tissue shows no (or very little) phase change. The amplitude 
information can be used to modify the phase correlation output to reject the noise, e.g. 
when the amplitude is small, the phase may be meaningless. The final output combines the 
gray scale amplitude information and color temporal phase information to display the B- 
5 mode with blood flow. This method is in contrast to the conventional color flow map, 
which relies on the Doppler effect of the flow, or other art that uses amplitude correlation 
for B-flow. The described phase application in B-mode, while not computing precise phase 
information at each sample position in each ultrasound frame, nevertheless produces useful 
flow information by taking advantage of the temporal distinctions between two successive 
10 frames. This is in contrast to other art which processes the phase information along the 
axis of the beam, as one is in the spatial domain and the other is in a fiiame-based temporal 
domain. 

Whereas many alterations and modifications of the present invention will no doubt 
become apparent to a person of ordinary skill in the art after having read the foregoing 
15 description, it is to be understood that the particular embodiments shown and described by 
way of illustration are in no way intended to be considered limiting. Therefore, reference 
to the details of the preferred embodiments are not intended to limit their scope, which is 
limited only by the scope of the claims set forth below. 
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CLAIMS 

What is claimed is: 

1 . A demodulator for use in an ultrasound infonnation processing system, said 
demodulator for receiving an input signal comprising digital samples, the digital samples 
5 being taken at a sampling frequency from an ultrasound transducer output, said 
demodulator for generating a demodulated output signal from said input signal, said 
demodulator comprising: 

a quadrature mixing device, said quadrature mixing device having a mixing 
frequency of one quarter of said sampling frequency, said quadrature mixing device 
10 comprising a half-band filter having a cutoff frequency of one-quarter of said 

sampling frequency, said quadrature mixing device for performing mirror- 
cancellation on said digital samples to produce intermediate signals; 
and 

a complex mixer for rotating said intermediate signals to a baseband 
15 frequency to produce an in-phase component and a quadrature phase component of 

said demodulated output signal; 

whereby an amplitude of said demodulated output signal corresponds to the vector 
magnitude of said in-phase component and said quadrature component. 

20 2. The demodulator of claim 1, said ultrasound transducer generating acoustic bursts 
at a carrier frequency, wherein said complex mixer has a mixing frequency equal to one 
quarter of the sampling frequency minus a first frequency function, whereby said 
demodulated output signal contains frequency components substantially near a frequency 
equal to said carrier frequency minus said first frequency function. 

25 

3. The demodulator of claim 2, said half-band filter of said quadrature mixing device 
being a first half-band filter, said quadrature mixing device comprising: 
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a first mixer for multiplying said digital samples with a first sinusoid, said 
first mixer having an output coupled to an input of said first half-band filter; 

a second mixer for multiplying said digital samples with a second sinusoid 
in quadrature phase with said first sinusoid; 
5 and 

a second half-band filter, said second half-band filter having an input 
coupled to an output of said second mixer; 

wherein said intermediate signals comprise an in-phase intermediate signal at the 
output of said first half-band filter and a quadrature-phase intermediate signal at the output 
10 of said second half-band fiher. 



4. The demodulator of claim 3, wherein said first mixer, said second mixer, said first 
half-band filter, and said second half-band filter are implemented on a single hardware 
chip, and wherein said first and second half-band filters are implemented in said hardware 

15 chip by finite impulse response filters having at least 40 taps. 

5. The demodulator of claim 4, wherein said sampling frequency is greater than 40 
MHz. 

20 6. The demodulator of claim 2 wherein, an a Doppler processing mode, said first 
firequency function is fixed at said carrier frequency. 

7. The demodulator of claim 6 wherein, in a B-mode processing mode, said first 
fi-equency function is a swept frequency function between a minimum fi-equency value and 

25 a maximum firequency value. 

8. A demodulator for use in an ultrasound information processing system, said 
ultrasound system comprising an ultrasound transducer generating acoustic bursts at a 
carrier fi-equency, the ultrasound transducer receiving acoustic echo signals and generating 

30 an input signal therefi-om, the ultrasound information processing system digitizing the 
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input signal at a sampling frequency to generate digital samples, said demodulator 
comprising: 

a quadrature mixer adapted to multiply said digital samples with 
quadrature-phase sinusoids at one quarter of said sampling frequency and filter the 
5 results with half-band filters each having a cutoff frequency of one-quarter of said 

sampling frequency to produce mirror-canceled signals; 
and 

a complex mixer for rotating said mirror-canceled signals to baseband to 
produce complex mixer output signals, whereby amplitude and phase information 
10 of a demodulated signal can be derived from said complex mixer output signals. 

?. The demodulator of claim 8, fiirther comprising: 

a secondary local oscillator for producing mixing signals that are mixed 
with the mirror-canceled signals by the complex mixer, said mixing signals being 
15 sinusoids at a frequency eqxial to the carrier frequency minus a first frequency 

function; 
and 

for each complex mixer output signal, a programmable low-pass filter; 

wherein, in a harmonic imaging mode, said first frequency function is set 
20 substantially near a frequency equal to twice the carrier frequency; 

and wherein, in said harmonic image mode, said programmable low-pass 
filter is adapted to have a notch at the carrier frequency for producing a 
programmable low-pass filter output having suppressed first harmonic information 
and passed second harmonic information &om said input signal. 

25 

1 0. The demodulator of claim 9, further comprising, for each programmable low-pass 
filter, an integrate-and-dump matched filter coupled to the output thereof for enhancing 
signal to noise perfonnance of said demodulator. 
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1 1 . The demodulator of claim 1 0, said demodulator also having a Doppler processing 
mode, wherein in said Doppler processing mode said first frequency function of said 
secondary local oscillator is fixed at said carrier fi-equency. 

5 12. The demodulator of claim 1 1 , said demodulator also having a B-mode processing 
mode, wherein in said B-mode processing mode said first fi-equency function is a swept 
fi^quency function between a minimum fi*equency value and a maximum frequency value. 

13. The demodulator of claim 12, wherein said quadrature mixer is implemented on a 
10 first hardware chip having an at-least 40 tap finite impulse response (FIR) half-band filter, 

and wherein said programmable low-pass filter is implemented on a second hardware chip 
having substantially fewer taps than said first hardware chip, whereby said demodulator 
has reduced implementation cost by using fixed half-band filters for mirror canceling while 
using programmable digital filters in locations that do require a high number of FIR filter 
15 taps. 

14. The demodulator of claim 13, wherein said sampling frequency is greater than 40 
MHz. 

20 15. A method for demodulating an ultrasound signal, said ultrasound signal having a 
carrier fi-equency, comprising the steps of: 

sampling said ultrasound signal at a sampling fi-equency to produce digital samples; 
mixing said digital samples with quadrature-phase signals at frequency equal to 
one-quarter of the sampling frequency; 
25 filtering results of said mixing step with half-band filters each having a cutoff 

frequency at one-quarter of the sampling frequency to produce mirror-canceled signals; 
and 

rotating said mirror-canceled signals to baseband using a complex mixer. 

30 
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16. The method of claim 15, wherein said sampling frequency is at least four times 
greater than said carrier frequency. 

17. The method of claim 1 6, further comprising the step of, prior to the step of 
5 sampling, filtering said ultrasound signal with an anti-aliasing filter having a cutoff 

frequency not greater than half the sampling frequency. 

1 8. The method of claim 1 7, said step of rotating comprising the step of mixing said 
mirror-canceled signals with quadrature-phase sinusoids having a frequency equal to one- 

10 fourth of said carrier frequency minus a first frequency function, said first frequency 
function being equal to said carrier frequency in a Doppler mode, and said first frequency 
function being a swept frequency function a B-mode. 

1 9. The method of claim 1 8, said first frequency function being substantially near twice 
IS the carrier frequency in a harmonic imaging mode. 

20. The method of claim 19, further comprising the step of, in said harmonic imaging 
mode, low-pass-filtering at least one of said rotated signals with a low-pass filter having a 
notch at the carrier frequency for producing an output having suppressed first harmonic 

20 information and passed second harmonic information from said input signal. 

21. An imaging method comprising: 

sending ultrasoimd energy into a body, receiving from the body ultrasound 
energy related to interactions between the body and the ultrasound energy sent 
25 therein, and deriving an input signal related to the received ultrasound energy; 

converting the input signal into digital samples related to a sampling 
frequency; 

processing the digital samples through a digital quadrature mixer using first 
mixing signals at a frequency of approximating a quarter of the sampling frequency 
30 to produce quadrature mixer outputs; 



/ 
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10 



filtering the quadrature mixer outputs through at least one digital filter 
having a passband approximately half the sampling frequency to produce 
intermediate digital signals; 

processing the intermediate digital signals through a complex mixer to 
produce complex mixer outputs; 

processing the complex mixer outputs to derive estimates of information 
related to at least one of an envelope of the input signal and a phase of the input 
signal; and 

using the estimates to generate and display an image related thereto. 

22; A method as in claim 21 in which the filtering comprises filtering the mixer outputs 
through at least one digital half-band filter. 



23. A method as in claim 22 in which the processing of the complex mixer outputs 

15 comprises processing the complex mixer output through at least one low pass filter 

to produce low pass filter outputs, 

24. A method as in claim 22 including processing the low pass filter outputs through at 
least one integrate-and-dump matched filter. 

20 

25. A method as in claim 21 in which the processing of the intermediate signals 
through the complex mixer comprises mixing the intermediate signals with second 
mixing signals at a frequency approximating the algebraic sum of quarter the 
sampling frequency and a frequency fimction. 

25 

26. A method as in claim 25 in which the imaging comprises a Doppler mode, and the 
frequency fimction is fixed. 

27. A method as in claim 25 in which the imaging comprises a B-mode and the 
30 frequency fimction varies in a manner related to changes in an instantaneous 

frequency of the input signal. 
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28. An ultrasound imaging system comprising: 

a transducer sending ultrasound energy into a body and receiving from the 
body ultrasound energy related to interactions between the body and the ultrasound 
5 energy sent therein, said transducer producing an input signal related to the 

ultrasound energy received thereby; 

a converter receiving the input signal produced by the transducer and 
deriving therefrom digital samples related to a sampling frequency; 

a source of first mixing signals at a frequency approximating a quarter of 
10 the sampling frequency; 

a digital quadrature mixer receiving the first mixing signals and the digital 
samples and in response producing quadrature mixer outputs related thereto; 

at least one digital filter receiving and filtering the quadrature mixer outputs 
to provide a passband approximately half the sampling frequency and thereby 
15 produce intermediate digital signals; 

a source of second mixing signals; 

a digital complex mixer receiving the second mixing signals and the 
quadrature mixing outputs to produce complex mixer outputs; 

a detector receiving and processing the complex mixer outputs to derive 
20 estimates of information related to at least one of an envelope of the input signal 

and a phase of the input signal; and 

a display displaying information related to the estimates of at least one of 
the envelope and the phase. 

25 29. A system as in claim 28 in which the digital filter comprises at least one digital 
half-band filter. 

.30. A system as in claim 28 in which the complex mixer comprises at least one low 
pass filter produce low pass filter outputs, 

30 
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A system as in claim 30 including at least one integrate-and-durap matched filter 
processing the low pass filter outputs. 

A system as in claim 28 in which the source of the second mixing signals 
comprises a source of signals at a frequency approximating the algebraic sum of 
quarter the sampling frequency and a frequency fiinction. 

A system as in claim 32 in which the source of the second mixing signals produces 
signals at fixed frequency for system operation in a Doppler mode but at a 
frequency varying in a manner related to changes in an instantaneous frequency of 
the input signal for system operation in a B-mode. 



-26- 




1/13 



wo 01/54583 



PCT/DSOO/35723 




2/13 



wo 01/54583 PCT/USOO/35723 




CM 

CD 
it 



3/13 



wo 01/54583 



PCT/USOO/35723 




5/13 



wo 01/54583 



PCT/USOO/35723 




6/13 



wo 01/54583 



PCT/USOO/35723 



< 


sin 


o 

CO ,.>c 






< 








cos 


-J. 

i i. 




ej 



§ 



7/13. 




8/13 



wo 01/54583 




PCT/USOO/35723 



5- 




9/13 



wo 01/54583 



PCT/USOO/35723 




10/13 



wo 01/54583 



PCT/USOO/35723 




11/13 



wo 01/54583 



PCT/USOO/35723 





/ 







u 


























i 




' — i 


: — ' 







CO 


Co X 


4 










< 


COS 


o it 
-J 




CD 



>8 



12/13 



wo 01/54583 



PCT/USOO/35723 




13/13 



INTERNATIONAL SEARCH REPORT 



International application No. 
PCT/US00/3S723 



A. CLASSIFICATION OF SUBJECT MATTER 

IPC(7) :A61B 08/00 
US CL :600/443 

According to International Patent Classification (IPC) or to both national classification and IPC 



B. FIELDS SEARCHED 



Minimum documentation searched (classification system followed by classification symbols) 
U.S. : 600/443 



Documentation searched other than minimum documentation to the extent that such documents are included in the fields searched 
600/447 



Electronic data base consulted during the international search (name of data base and, where practicable, search terms used) 



DOCUMENTS CONSIDERED TO BE RELEVANT 



Category* 



Citation of document, with indication, where appropriate, of the relevant passages 



Relevant to claim No. 



US 5,437,281 A (UN et aL) 1 August 1995. See entire document. 

US 5,445.156 A (DAFT et al.) 29 August 1995. See entire 
document. 



1-33 
1-20 



\ I Further documents are listed in the continuation of Box C. See patent family annex. 



Special calcgories of cited docunienu; 

document defining the general state of the art which ii not eoiuidered 
to be of panicular relevance 

earlier document published on or after the intBrnationBl fiUng dale 

document which may throw doubU on priority claiin(s) or which ii 
cited to estBblish the publication date of another citation or other 
special reason (as specified) 

document referring to an oral disclosure, use. exhibition or otho- 



documcnt published prior to the intemattonal filing dale but later than 
(he priority dale claimed 



*T* later document published after the international niing dale or priority 

date and not in conflict with the application but cited to understand 
the principle or theOl> undeiiying the invention 

*X' document of paiticular relevance; the claimed invention cannot be 

eoiuidered novel or cannot be considered to involve an inventive step 
when the document is taken alone 

'Y* doeumcnl of particular relevance; the claimed invention cannot be 

considered to involve an inventive step when the document is 
combined with one or more other such documents, such combination 
being obvious to a person skilled in (he art 

*&' document member of the same patent family 



Date of the actual completion of the international search 



12 FEBRUARY 2001 



Name and mailing address of the ISA/US 
Commissioner of Patents and Trademarks 
Box PCT 

Washington, D.C 30231 
Facsimile No. (703) 305-3230 



Date of mailing of the international search report 




3 3 MAR 2001 



A\uHbpz^of|^cer/'^N 
MAULIN PATEC ' 
Telephone No. (703) 305-6933 




Form PCT/ISA/210 (second sheet) (July 1998)* 



wo 01/54583 



PCT/USOO/35723 



mirror-cancelled signals. The demodulator further includes a complex mixer for rotating 
the mirror-canceled signals to baseband, whereby amplitude and phase information of the 
desired demodulated signal can be derived from said complex mixer output signals. The 
complex mixer is provided with mixing signals from a secondary local oscillator, the 
5 mixing signals being sinusoids at a frequency equal to the carrier frequency minus a first 
frequency function. For Doppler mode imaging, the first frequency function is fixed at the 
earlier frequency. For B-mode-imaging, the first frequency function may be a swept 
frequency function. 

In accordance with a preferred embodiment, a programmable low-pass filter is 
10 provided at each output of the complex mixer, and an integrate-and-dump matched filter is 
provided at the output of each programmable low-pass filter. Advantcigeously, a harmonic 
imaging mode may be achieved by setting the first frequency function of the secondary 
local oscillator substantially near a frequency equal to twice the carrier frequency, and by 
adjusting the progranmiable low-pass filter to have a notch at the carrier frequency. Using 
15 these settings, the output of each programmable low-pass filter comprises a signal 
containing second harmonic information and suppressing first harmonic information. 

BWEF DESCRIPTION OF THE DRAWING 

FIG. lA illustratesanexampleof an ultrasound transducer input signal; 
20 FIG. IB illustrates an example of an ultrasound transducer output signal; 

FIGS. IC illustrates frequency spectra of portions of the output signal of an 
ultrasound transducer corresponding to differing target depths; 

FIG. 2 is a block diagram of an ultrasound information processing system in 
accordance with prior art; 
25 FIG. 3 is a block diagram of a demodulator used in the ultrasound information 

system of FIG. 2; 

FIG. 4 comprises parts 4-1 and 4-2 and illustrates frequency spectra corresponding 
to the operation of the demodulator of FIG. 3; 

FIG. 5 is a block diagram of a demodulator for use in an ultrasound information 
30 processing system in accordance with a preferred embodiment; 



